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ABSTRACT 

Although a source star is fainter than the detection limit imposed by 
crowding, it is still possible to detect an event if the star is located in the seeing 
disk of a bright star is and gravitationally amplified: amplification bias. Using 
a well-constrained luminosity function, I show that ~ 40% of events detected 
toward the Galactic bulge are affected by amplification bias and the optical 
depth might be overestimated by a factor ~ 1.7. In addition, I show that if one 
takes amplification bias into consideration, the observed time scale distribution 
matches significantly better, especially in the short time-scale region, with 
the distribution expected from a mass-spectrum model in which lenses are 
composed of the known stellar population plus an additional population of 
brown dwarfs than it is without the effect of the amplification bias. 
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1. Introduction 



Many candidates of Massive Astro-physical Compact Halo Objects (MACHOs) are 
being detected by the MACHO (Alcock et al. 1993, 1996a), EROS (Aubourg et al. 1993, 
1995), OGLE (Udalski et al. 1992, 1994), and DUO (Alard 1996) groups. To maximize 
the event rate, the searches are being carried out toward very dense star fields, e.g., the 
Galactic bulge and Magellanic Clouds, in which the detection limit is set by crowding. 
However, it is possible to detect a lensing event with a source star that is below the 
detection limit, and thus unresolved, provided that the star is located in the seeing disk of 
a bright star and is gravitationally amplified. This effect is known as "amplification bias" 
(Blanford k, Narayan 1992; Narayan & Wallington 1994). I will call events of this type and 
the corresponding source stars below the detection limit "faint events" and "faint stars", 
respectively. Current lensing experiments are adopting an observational strategy in which 
they construct templates of resolved stars to allow fast photometric comparison of source 
star luminosities. It may then appear as if faint events cannot be detected because they 
are not registered in the template (Bouquet 1993). However, in the very dense field in 
which nearly all stars are blended, the amplified flux of a faint star will increase the flux 
of a nearby bright star, and thus an event seemingly with the bright star as a source star 
can be detected. I will call the stars registered in the template image "bright stars" as 
opposed to "faint stars" . Therefore, the current lensing searches suffer from amplification 
bias (Nemiroff 1994). 

In this paper, using a well-constrained luminosity function (LF) I show that ~ 40% 
of events detected toward the Galactic bulge are affected by amplification bias and that 
the optical depth might be overestimated by a factor ~ 1.7. In addition, I show that by 
taking amplification bias into consideration, the observed time scale distribution matches 
significantly better, especially in the short time-scale region, with the distribution expected 
from a mass spectrum model in which the lenses are composed of the known stellar 
population plus an additional population of brown dwarfs than it does without taking 
account of amplification bias. 
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2. Faint Lensing Events 



For the detection of a faint event with a faint star flux L < Ldl, the amplification 
should satisfy the condition 

F_ L b + A min L _3_ 

F L b + L ^5' 1 j 

resulting in the minimum required amplification of 

A m ,„= (3/V5)( Y L) -^ (2,2) 

Here L& and L_d_l are the fluxes of the blended bright star and the detection limit. The 
factor 3/y/b in above equations is included so that an event can be detected when it is 
amplified more than A = 3/y/b of the combined flux of the bright and faint stars. The 
amplification is related to the lensing parameters by 

u 2 + 2 2 a2 (t-t Q ) 2 

A{u) = u(u 2 + ^ U= ^^T' ( 3) 

where u is the lens-source separation in units of the Einstein ring radius r e , (5 is the impact 
parameter, to is the time of maximum amplification. The Einstein ring radius is related to 
the physical parameters of a lens by 

where Ml is the mass of the lens, and D Q \, D\ s , and D os are the distances between the 
observer, lens, and source. The Einstein time scale is related to r e by t e = r e /v, where v is 
the transverse speed of the lens relative to the observer-source line of sight. The maximum 
allowed impact parameter for detection is related to the minimum required amplification 
by 
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For faint events, one can detect only the portion of a light curve above the threshold Lf, 
and the event therefore mimics that of a shorter time-scale event. Then what is measured 
is not t e but the effective time scale t e flf ; which is the half of the duration of an event above 
the threshold. The ratio between these two time scales is 

^) = ^=[/?Lx-/3 2 ] 1/2 <l-0. (2.6) 
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3. The Fraction of Amplification-biased Events 

Due to amplification bias, each individual bright star works effectively as multiple 
source stars. Let B the effective number of source stars per single bright star: amplification 
bias factor. This factor is computed by 

B(L b ) = 1 + dL$(L)(3 max (L b , L){rj), (3.1) 
Jo 

where ®(L) is the LF of Galactic bulge stars normalized to the area of the seeing disk 
around each bright star. The average seeing of the current experiments is ~ 2". However, 
when a faint source is located at the edge of the seeing disk of the bright star, one can 
isolate its magnified image from the bright star. I therefore set the undistinguishable 
separation between images at AO = 1".5, i.e., §{L) is normalized for stars in the area 
7r(l".5) 2 . In equation (3.1) the value /9 max is included because only events with j3 < /3 max 
can be detected. In addition, the factor (77) = J(f max r)((3) d(3/f3 maiX enters because the 
detection efficiency, e, decreases as the time scale decreases. For the moment I assume 
that the efficiency is linearly proportional to the time scale. The adopted values of A9 and 
the functional form e(t c ) is subject to some uncertainties, so I will discuss other cases in § 
4. Then the total effective number of source stars is obtained by integrating B(L b ) for all 
bright source stars weighted by the LF; 

[•OO 

B tot = / dL b $(L b )B(L b ) 

poo roo rLfj 

= dL b <f>(L b ) + dL b <$>{L b ) dL<f>(L)(3 max (L b ,L)(ri). (3.2) 
JLol jLdl Jo 

The first and second terms in equation (3.2) are the number of bright stars in the 

seeing disk and the additional faint stars that effectively work as source stars due to the 

amplification bias effect, respectively. 

The model LF is constructed as follows. For stars brighter than the de-reddened 
/-band mag of Iq = 18.2, I adopt the LF determined by J. Frogel (1996, private 
communication), and in the range 18.2 < Iq < 22.4, I use the LF determined by Light 
et al. (1996). To extend the LF beyond even this limit, I adopt the LF of stars in the 
solar neighborhood determined by Gould, Bahcall, & Flynn (1996). I address below the 
uncertainty caused by the difference in the bulge and disk stellar populations. I adopt a 
distance i?o = 8.0 kpc to the Galactic bulge stars. The model LF is presented in Figure 
1 in units of stars mag -1 arcmin -2 . Current experiments reach the detection limit when 
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the stellar number density of Galactic bulge fields arrives at ~ 10 6 stars deg 2 (C. Alcock 
1996, private communication). Based on the model LF this number density corresponds to 
I = 18.2 mag or Mj = 3.7. 

With the model LF and the corresponding detection limit, the amplification bias 
factor B is computed by equation (3.1) along with equations (2.2), (2.5), and (2.6). [] The 
resultant values of B are shown as a function of bright star mag in Figure 2. Because the 
LF increases as Lb decreases, combined with the fact that the most probable faint events 
have unamplified faint-star flux just below L b , the value B increases as Lb decreases. Since 
the event rate is directly proportional to the number of source stars, the average increase 
in event rate due to the amplification bias effect is determined by 

(B) = roc y (3-3) 

Jl dl dL h $>{L b ) 

resulting in (B) = 1.65. That is, on average each source star works effectively as ~ 1.7 
stars, and thus the event rate increases by the same factor. 

The distribution of faint source-star brightness for a given template star with a flux 
Lb is obtained by taking the derivative of B, i.e., 

dB{ ^ L) = S(L - L b ) + P max (L b , L)( V ML). (3.4) 

Note that dL5{L — Lb) = 1. The distributions dB(Lb,L) for bright stars with 
lb = 15, 16, 17, 18 mag are shown in the left panel of Figure 3. The total event 
distribution for all bright stars is computed by 

/■OO 

dB tot {L)= dL b <5>{L b )dB{L b ,L) 

J Lr>L 

/■oo 

= $(L b )dL b + / dL b <S>(L b ) [P max (L b , L)( V )<S>(L)dL] , (3.5) 

where the first and second terms represent the event rate distribution without and with the 
amplification bias effect, respectively. When the amplification bias effect is not included, 



An additional threshold for detection might be imposed by the flux from other stars located in the 
seeing disk, background flux. However, flux from stars much fainter than Lb will spread smoothly over the 
seeing disk and it will be subtacted during image process. Some stars just below Lf, that are located in the 
seeing disk will cause the threshold to be higher. On the other hand, stars that are near but not inside the 
seeing disk will make the bright star appear fainter and so cause a lower threshold. To lowest order, these 
two effects cancel one another. 
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the distribution is just proportional to the LF, i.e., only the first term. The total event 
distribution is shown in the right panel of Figure 3, in which the contribution to the 
distribution by faint events are shaded. The faint events comprises ~ 40% of the total 
events and the unamplified mag of faint source stars extends up to Iq ~ 23. For very faint 
stars, i.e., Iq ^ 23, to be detected they should be highly amplified, and thus they are 
rare. The slim chance of high amplification becomes even slimmer because of the drop in 
efficiency. Therefore, the uncertainty in the model LF in the very faint region does not 
affect the determination of (£>). 



4. Effect of the Amplification Bias on the Time Scale Distribution 



Amplification bias has a major influence on the Einstein time scale distribution f(t e ) 
because one measures t e g instead of t e . How and how much is the time scale distribution 
affected? To answer this question I construct time scale distributions for events toward 
the Galactic bulge with and without including the effect of amplification bias using a 
reasonable mass spectrum model of lenses. For the construction of f{t e ), it is required 
to model the mass and velocity distributions. For the Galactic bulge mass distribution, I 
adopt a "revised COBE" model that is based on the triaxial COBE model (Dwek et al. 
1995) except for the central part of the bulge. In the inner ~ 600 pc of the bulge, I adopt 
the centrally concentrated axisymmetric Kent model (1992) since the COBE model does 
not match well in this region. The model is provided in terms of the light density, v. For 
the disk, I adopt a Gould et al. (1996) model which has of the form 



p(R,z) = p 



r ech2 (£> + o" exp (-xr 



'2 



where h\ = 323 pc and hi = 660 pc, and po = 0.436 M Q pc~ 3 . Both the disk 
and bulge MACHO transverse speed is modeled by a Gaussian. In the model, 
the velocity distributions of disk MACHOs have means and standard deviations of 
(v y ,a y ) = (220,30) km s _1 and (v z ,a z ) = (0,20) km s _1 . The projected components of 
the Galactic bulge velocity dispersion are computed from the tensor virial theorem and 
results in (v y ,a y ) = (0,93.0) km s" 1 and (v z ,a z ) = (0,78.6) km s _1 . Here the projected 
coordinates (y, z) are set so that the axes are respectively parallel and normal to the 
Galactic plane. For more details of the mass and velocity distributions, see Han & Gould 
(1995, 1996). 



-7- 



There may be dark lenses as well as lenses from known stellar populations. The mass 
spectrum, /(Ml), of the stellar population is constructed by using the mass-luminosity 
relation provided by equation (5) of Henry & McCarthy (1993). The white dwarf 
component in the mass range 0.5 Mq < Ml < 0.7 Mq is included by normalizing its mass 
spectrum so that there are ~ 10 times more white dwarfs than the number of turnoff 
plus giant stars. Finally, brown dwarfs in the mass range 0.07 Mq < Ml < 0.09 M & 
are included in the mass spectrum making up the rest of the total bulge mass of 
Afbuige = 2.1 x 10 10 M o , which is adopted from Zhao, Spergel, & Rich (1995). Then the 
total mass of each lens population i in the Galactic bulge is determined by 

M pop , = M BW , (^) , (4.2) 

where L b ui g e = J hulge dxdydz v(x,y,z) = 1.8 x 10 10 L© and L B w = Je BW dD ol v(D ol ) = 
2412 L© pc~ 2 is the total amount of light in the bulge and the integrated light seen 
through a unit area (pc 2 ) of the Baade's Window (BW) (Kent 1992; Dwek et al. 1995). 
The integrated mass of each population is obtained from the mass spectrum model, i.e., 
Mbw,! = / dM L fi(M L )M L . With the mass spectrum model I find that M BW ,i = 1579, 339, 
and 897 Mq pc~ 2 , and the resulting total masses in the bulge are M poPi j = 1.18 x 10 10 , 
0.25 x 10 10 , and 0.67 x 10 10 Mq for the stellar, white dwarf, and brown dwarf populations, 
respectively. 

With these models, the event rate distribution of bulge-bulge self-lensing and 
disk-bulge events for a single source star is computed by equations (3.3) and (3.4) of Han 
& Gould (1996). The lens masses are drawn from the model mass function. However, the 
COBE bulge model is given in terms of light density, u, and thus the conversion from v to 
n is required for proper normalization. This number-to-light density ratio is determined by 

n - = ^w. (4.3) 

Here Nbw = J2i ^BW,i is the total number of objects in the unit area of sky toward BW, in 
which the number of individual component is obtained by JVbw.s = / dML U^Ml) = 3441, 
566, and 11357 objects pc~ 2 for the stellar, white dwarf, and brown dwarf populations, 
respectively. Once again the total amount of light in the same area of sky is 
Lbw = 2412 Lq pc -2 [see below eq. (4.2)]. I find njv = 6.31. Note that although brown 
dwarfs comprise only 32% of the total mass, they account for 74% of the total number 
density. 

Once f(t c ) is obtained, the total event rate as a function of time scale not including the 
amplification bias effect, F w / Q (t e ), by monitoring all source stars is obtained by multiplying 
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the total number stars, A*, and the total amount of observation time, T, into f(t c ); 

( \ -^cvcnt j. 
E J^-y ( 4 - 4 ) 

where N event is the actually detected number of events. Note that iV* is the number of stars 
in the template, and thus only resolved stars. The MACHO group reported -/V event s = 39 
events in the first year bulge season, resulting in an optical depth of r = 2.4 x 10 -6 for all 
types of stars (Alcock et al. 1996b), which implies iV*T = 1.6 x 10 9 days. 

On the other hand, the construction of the total time scale distribution with 

amplification bias, Fwith^cff)) requires additional processing which is described below. For 

a single bright star, the factor rj is distributed as 

f°° f Lb /-/3max r , x i/2] r r°° l - 1 

g( V )= dL b <S>(L b ) dL<S>(L) d/3 6 r, - (p 2 max - (5 2 ) / dL b $>{L b ) . 

Jl dl Jo Jo L v 7 J Ul dl J 

(4.5) 

Once g(rj) is obtained, the effective time scale distribution of faint events Ff(t e g) is 
obtained from F w / (t c ) by 

p 1 /*oo 

F f (t eS ) = J o drjg(ri) J dt efl F w/o (t e )5(t e s - rjt e ). (4.6) 

Then, one finds the total (both bright and faint) event time scale distribution by 
-^with(^e) = Fw/o(te) + Ff(t e ). With the detection efficiency e(t e ) provided by Alcock et 
al. (1996b), the final time scale distributions with and without amplification bias are 
computed by r w/o (i c ) = e(t c )F w/o (t c ) and r with (t e ) = e(t c )F with (t e ), and they are shown 
in Figure 4. In the figure, the distributions are compared with the observed time scale 
distribution (shaded histogram) obtained by the MACHO group (Alcock et al. 1996b). 

There are two major changes in the time scale distribution by taking the amplification 
bias effect into account. First, the distribution r w j t h(i e ); in general, has a higher 
normalization relative to T w / Q (t e ) due to the increase in number of stars that are effectively 
monitored. The increase factor is (£>) = / dt c r w i t h(ie)/ / dt c r w/ / (t c ) = 1.65, which matches 
well with the value determined in § 3 using the approximation e(i c ) oc t c . Second, r w ; t h(£ e ) 
is shifted toward shorter time scale compared to T w / Q (t c ) due to the additional contribution 
by faint events which mainly have short time scales. Therefore, a significant fraction of 
short events (~ 10 days), which could not be explained by known lens populations, might 
be caused by the amplification bias effect. 

However, there are some uncertainties in determining the contribution to the event 
rate by faint events. One of the uncertainties comes from the maximum size of separation 
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at which one can isolate the image of the amplified faint star from that of a bright star. I 
compute the increase factors in event rate for A9 = 1".0 and 2".0 and find that (B) = 1.35, 
and 2.38, respectively. Additional uncertainty comes from the blending of bright stars 
with faint stars. Due to the blending, some fraction of bright events would fail to be 
detected and the measured time scale would be shorter, resulting lower normalization and 
additional shift toward shorter time scale in T(t e ). However, because of the dominance of 
the bright star flux, this effect would also be small. 



5. The Effects of Amplification Bias on Optical Depth Determination 

When the amplification bias is not taken into consideration the optical depth might 
be overestimated by a factor (B) ~ 1.7 because events are detected by monitoring 
A* jC fr = (B)N*, while r is determined with A* instead of A* iC g. On the other hand, 
the effects of amplification bias on t c and e do not propagate to the optical depth 
determination. This is because by measuring t e g instead of t e , the time scale decreases by 
a factor ~ r] and for the same reason the detection efficiency decreases by a similar factor, 
i.e., the same factors cancel each other out and thus there is no net effect. In addition, 
most additional events due to amplification bias are expected to have short t e s in which 
the efficiency is well approximated as linear. 



6. Conclusion 

I have shown that amplification bias may have significant effects on the time scale 
distribution and the optical depth determination, and thus the correction of the bias is 
very important. The true distribution T w / Q (t e ) can be recovered statistically with the 
known weight factor dBtot(L) in the reverse way that r w ; t h is obtained from r w / . In 
addition, for some long t e g events for which very detailed light curves can be constructed, 
it will be possible to recover the individual true t c by the fitting the light curves using 
an additional parameter, the unmagnified flux. Another way to detect biased events is 
finding the shift of the centroid of source stars caused by the faint star amplification. This 
shift was actually detected by Alard (1996) despite the moderate quality of the image: 
photographic plate. An inspiring development in the lensing experiments is that the time 
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rcsolution of observations is improving rapidly. For example, the alert system allows 
intensive network observations of candidate events (PLANET, Albrow et al. 1996; GMAN, 
Pratt 1996). The EROS group (Aubourg et al. 1995) has carried out a lensing experiment 
toward the Magellanic Clouds with a monitoring frequency of up to 46 times per night. 
In this way, it will be possible to obtain detailed light curves, and so t e for a significant 
fraction of individual events. 

I acknowledge precious discussions with A. Gould. This work was supported by the 
grant AST 94-20746 from the NSF. 
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Figure 1: The model luminosity function in de-reddened J-band. 
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Figure 2: The amplification bias factors as a function of bright source star mag. Bacause the LF of bright 
stars increases, combined with the fact that the most probable faint events have unamplified faint-star flux 
just below Lb, the value B increases as Lb decreases. 
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Figure 3: The distributions of faint source-star brightness for bright star mag of I b — 15, 16, 17, and 18. 
Also shown are the total (both bright and faint) event distributions as functions of source star mag. The 
contribution by faint events is shaded. Under the approximation that e « t e , faint events comprise 40% of 
total events. 
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Figure 4: The time scale distributions for all stars with and without including amplification bias effect, 
and they are compared with the observed distribution obtained by the MACHO group (shaded histogram). 



